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I. INTRODUCTION 

Lattice calculations of the kaon mixing parameter Bk 
are now very precise, with results using several types 
of fermion showing reasonable consistency. These re- 
sults play an important role in constraining the param- 
eters of the Standard Model (SM) [TJi]. Here we con- 
sider matrix elements of operators having Dirac struc- 
tures other than the "left-left" form that arises in the 
SM. These new operators generically appear in models of 
physics beyond the SM (BSM) when heavy particles (e.g. 
squarks and gluinos in supersymmetric models) are inte- 
grated out (see, e.g., Refs. [SHE]). The resulting AS = 2 
four-fermion operators give rise to additional contribu- 
tions to both the CP-conserving and CP-violating kaon 
mixing matrix elements, both of which are strongly con- 
strained by experimental results. Thus, if one can cal- 
culate the corresponding hadronic matrix elements, one 
can place significant constraints on BSM physics (see, 
e.g., Refs. USEE]). 

As part of the Staggered Weak Matrix Element 
(SWME) collaboration, we are undertaking a numeri- 
cal calculation of these matrix elements using improved 
staggered fermions — specifically, HYP-smeared [3] va- 
lence quarks on lattice configurations generated by the 
MILC collaboration with Nf = 2 + 1 flavors of asq- 
tad sea quarks [TO]- To extrapolate or interpolate the 
results to the physical d and s quark masses and the 
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continuum limit without introducing model dependence, 
it is advantageous to use functional forms incorporat- 
ing the constraints of chiral symmetry. These forms can 
be calculated using staggered chiral perturbation theory 
(SChPT) [TThT5] . and the present paper provides the re- 
sults at next-to-leading order (NLO). One also needs to 
match the lattice operators onto those regularized in a 
continuum scheme, and the necessary matching factors 
were previously calculated to one-loop order in pertur- 
bation theory [14]. 

The corresponding analysis for the Bk operator has 
been carried out in Ref. |15j . and turns out to be quite 
challenging. The difficulty arises because the left-left 
structure of the four-fermion operator leads to suppres- 
sion of the K — K matrix element in the chiral limit. 
However, many of the operators which arise from mix- 
ing due to discretization errors and truncated perturba- 
tive matching do not have chirally suppressed matrix el- 
ements. This leads to a plethora of unknown low-energy 
constants (LECs) entering at NLO in SU(3) SChPT. The 
situation is much simpler, however, in SU(2) SChPT, 
where there is only one additional LEC at NLO com- 
pared to continuum ChPT. 

Similarly, the results for the BSM four-fermion opera- 
tors are much simpler than for the Bk operator because 
none of the BSM operators have chirally suppressed ma- 
trix elements. In fact, the situation for both SU(3) and 
SU(2) SChPT for the BSM operators is the same as that 
for Bk in SU(2) SChPT. As long as one considers ap- 
propriate ratios ("B factors"), there is only a single ad- 
ditional LEC compared to the continuum ChPT expres- 
sions. This new LEC is induced by matching and dis- 
cretization errors. The simplicity of the SChPT result 
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should make the extrapolation of the new matrix ele- 
ments straightforward. 

One-loop results in continuum ChPT for the BSM op- 
erators have been given in Ref. [16 . These provide an 
important check on our results. As a spin-off from our 
calculation, we provide the partially quenched generaliza- 
tion of the results of Ref. pi] for both SU(3) and SU(2) 
ChPT. 

This paper is organized as follows. In Sec. [IT] we list the 
BSM operators and describe how they are mapped into 
the partially quenched (PQ) lattice theory. In Sec. Ill 
following a brief review of SChPT, we explain the map- 
ping of the lattice operators into the chiral effective the- 
ory. Section |TV] presents the NLO calculation and results, 
first for SU(3) PQSChPT and then for SU(2) PQSChPT. 
We also give the continuum SU(3) and SU(2) results in 
PQChPT. We close in Sec.fVlwith some brief conclusions. 



II. OPERATORS AND B-FACTORS 

A. Continuum operators 

We use the so-called SUSY basis for the BSM opera- 
tors pi] 



2 = s a (l- Tj )d a s b (l- l5 )d b , 

3 = s a (l- l5 )d»s»(l- l5 )d a , 
A = s a (l- l5 )d a s\l+ 75 )d b , 
5 = s a (l- l5 )d b s b (l+ l5 )d a , 



(1) 

(2) 
(3) 
(4) 



where a and b are color indices. Together with the Bk 
operator 



1 =s a 1 ^l- l5 )d a s b 1 ^l- l5 )d 



(5) 



they form a complete set of AS = 2 four-fermion opera- 
tors. In the following we will concentrate entirely on the 
BSM operators, and the index j will always run over the 
values 2 — 5. 

These operators must be renormalized, and the Wil- 
son coefficients which multiply them are usually calcu- 
lated in a canonical choice of continuum scheme (such 
as MS with naive dimensional regularization of 75) at a 
canonical scale (such as 2 GeV). Consequently, we are 
ultimately interested in the matrix elements of the oper- 
ators defined in such a scheme. On the lattice, however, 
one inevitably starts with bare lattice operators, or with 
operators defined in a regularization independent scheme 
such as RI-MOM, and one must match the operators to 
those in the canonical scheme. In our ongoing numerical 
work we presently use one-loop perturbative matching to 
bare lattice operators, using the results from Ref. [H] , 
Thus, in the following we assume that we are using those 
linear combinations of lattice operators which match at 
one-loop order to the operators defined in the canonical 
continuum scheme. 



It is convenient and conventional to package our igno- 
rance of the matrix elements of O2-5 into B-parameters. 
For the BSM operators these are defined as follows [T5] : 



_ {K \O 3 ^)\K ) 

N j (K \s^ 5 d a (fJ,)\0)(0\s'> l5 d b (^\K )' 



(N 2 ,N 3 ,N 4 ,N 5 ) = (5/3,-1/3,-2,-2/3), 



(0) 



(7) 



where [i is the renormalization scale. The denomina- 
tors in these ratios are obtained using the vacuum sat- 
uration approximation, including the contribution from 
the Fierz-rearranged form of the operators, but dropping 
contributions suppressed in the chiral limit. The product 
of matrix elements in the denominators can be written 



(ir |s7 5 ^)|0><0|s 75 <M)l#o) = - 



m d (n) + m s (/x) 



(8) 

which makes explicit that both the numerator and de- 
nominator in the ratios defining the Bj depend on the 
renormalization scale, and should be defined in a com- 
mon scheme. 

There are several advantages to using the ratios Bj 
rather than directly calculating the matrix elements 
(K a \Oj(fi)\K ). First, in a lattice calculation, forming 
a dimensionless ratio reduces statistical and systematic 
errors — particularly those due to the uncertainty in the 
lattice spacing and in the matching factors. Second, the 
SWME lattice calculation uses wall sources, following the 
same methodology as for Bk [H], and the overlap fac- 
tor between these sources and the kaon states cancels in 
the ratio. Third, as we will see below, the SChPT ex- 
pression for the ratio is simpler (involving fewer LECs) 
than for the matrix elements. There is, however, also 
a potential disadvantage, as stressed in Ref. [20]. To 
convert from the Bj to the corresponding matrix ele- 
ments, one must multiply by the denominator, which, as 
shown by Eq. ([8]), depends on nid and m s . These quark 
masses are not directly measurable physical quantities, 
and must be obtained from lattice calculations. In the 
late 1990's, when Ref. [2D] was written, there were large 
systematic errors in determinations of light quark masses; 
the source of largest error was the quenched approxima- 
tion. The present situation is markedly improved, with 
quark masses known to about 1% accuracy [TJ[2]. Thus, 
there is no longer a phenomcnological reason not to use 
the Bj. 



B. Lattice operators in the continuum limit 

To calculate the Bj with staggered fermions, we must 
account for the additional taste degree of freedom. Each 
quark flavor enlarges to a quartet with four tastes. In this 
subsection, we first consider the staggered theory in the 
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continuum limit, in which the taste symmetry is exactQ 
Taste-breaking corrections will be discussed in the next 
subsection. 

The additional tastes occur both in the operators Oj 
and in the external kaons. We choose the latter to have 
taste P, i.e. to be created by operators with spin-taste 
75 ® £5- This kaon is the pseudo-Goldstone boson (PGB) 
associated with the spontaneous breaking of an axial 
U(1)a symmetry which holds exactly on the lattice in 
the massless limit. It follows that its correlation func- 
tions satisfy Ward-Takahashi identities which are ana- 
logues of those in the continuum [21] . This in turn leads 
to simplifications in the SChPT expressions for its matrix 
elements. The taste P kaon is also the simplest choice for 
numerical calculations, since it is the lightest kaon state. 

Turning now to the operators, we face the problem 
that Fierz transformations of the continuum Oj are no 
longer matched by those of the lattice operators once we 
introduce taste. Following Refs. [E1I22], we resolve this 
problem by introducing two types each of valence s and 
d quarks. We label these Si and S 2 (or D\ and D 2 ) using 
uppercase letters to denote fields which include the taste 
degree of freedom. Then the operators in the continuum 
staggered theory take the form [T3] 



/5-iCont /7-)Cont' , //-jCont 

U 2 — U S,II + U P,II 



, 1 



\ySJ + U P,/ ~ U T,I 



(9) 



U 3 - U S,I + U P,I - 2 \ U SJI + U P,II - U T,II J> 

(10) 

OCont' //-jCont' //-jCont' ^ / /oCont' /oCont'\ /i i \ 

4 — U S,II ~ U P,II — 2 \ U V,I ~ U A,I Ji [ iL ) 

OCont' />nCont' /nCont' /VoCont' /r\Cont'\ /in\ 
5 — U S,I ~ U P,I ~2 V v - n ~ U AJI )■ \ iZ ) 



Here the subscripts indicate firstly the "spin" of the four- 
fermion operator and secondly the manner in which the 
color indices are contractedrl The prime in the super- 
script Cont' is a reminder that this is the continuum the- 
ory in which the number of valence s and d quarks have 
been doubled. 

The "two-color-loop" operators (denoted by subscripts 



1 We are assuming here that the rooting procedure used to remove 
the additional tastes from the quark-gluon sea defines a theory 
with the correct continuum limit. 

2 The notation is slightly changed from that in Ref. 14 so as to 
conform to the more convenient notation of Ref. |15| . 



"IF) are 



O 



Cont' 
P,II 

TCont 



s?(75®6)£?^(75®&)r>2, 



(13) 
(14) 



O^fj = ^ St{ %llv ® S b 2 { llllv ® &)Dl, (15) 



O 



Cont' _ 
V.I I 



(16) 



O 



Cont' _ 
A, 1 1 — 



: ® S b 2 { llll5 ® (17) 



and are so named because, when contracted with exter- 
nal color-singlet kaon fields, there are two loops of color 
indices. The corresponding "one-color-loop" operators 
differ only in their color indices, as exemplified by 



^Cont' = ga {1 g, C5)D , ^ ^ )D a m 



(18) 



The matrix elements of the operators Gont are to be 
taken between a taste-P kaon of type 2, Kp 2 , created by 
the operator 1)2(75 ® £,5)82, and an antikaon of type 1, 
K PXl destroyed by ^1(75®^)^. In this way we force the 
component operators of the £)Cont t have a single Wick 
contraction with the external kaon fields, and thus avoid 
the Fierz-transformed contractions which occur in the 
matrix elements of the original operators Oj . The latter 
contractions are then added back by hand, giving rise to 
the terms in parentheses in Eqs. (9p2). Note that since 
the external kaons have taste P, the bilinears composing 
the four-fermion operators in Eqs. ( 13p8 ) must also have 
this taste, since taste is a good symmetry in the staggered 
continuum theory. 

The correspondence between i?-parameters in QCD 
and those in the augmented staggered theory can now 
be given. In the continuum limit of the latter theory, we 
have 



B 



2(K° P1 \Of ont '\K° 



P2/ 



^(^1^(75 ®&)A|0)(0|5 2 (75®6)£>2|tf'£ 2 )' 

In essence, we have constructed lattice operators which 
have the same Wick contractions with the external fields 
as do the original operators Oj between physical kaons. 
The extra factor of 2 in the numerator [compared to 
Eq. ([6|] accounts for the fact that in the original the- 
ory each bilinear could be contracted with either external 
field, whereas here there is only one such contraction due 
to the presence of two types of S and D quarks. In the 
staggered theory one also must account for possible fac- 
tors of the number of tastes, N t — 4. Such factors cancel 
in the ratios Bj |15j . and so are not shown explicitly. 

At this stage, it is helpful to summarize the content 
of the augmented staggered theory that we have con- 
structed. This theory contains U, D and S sea quarks, 
as well as D%, D 2 , Si and S 2 valence quarks. Each of 
these fields represents four degenerate tastes. We allow 
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the masses of the sea and valence quarks to differ, since 
we make use of this freedom in our simulations. We call 
the sea quark masses m u , and m a , respectively, while 
we follow Ref. [T5] and denote m^i = = rn x and 
m s i — m s2 — rriy. Note that we choose both strange 
valence quarks to have the same mass, and similarly for 
the down valence quarks. Finally, we must add ghost 
quarks for each of the valence quarks, and apply the 
fourth- root prescription to the sea-quark determinant. 
In the continuum limit, rooting is equivalent to adding 3 
tastes of ghost quark for each sea-quark field. Including 
both flavor and taste in the counting, the resulting par- 
tially quenched (PQ) theory has 28 quarks and 25 ghost- 
quarks. This is the minimal field content required to rep- 
resent the desired operators when using rooted staggered 
fermions. 

Although the construction of this PQ theory has been 
motivated by our use of staggered lattice fermions, one 
can also consider it as a purely continuum theory with no 
reference to the lattice. The result ( 19 1 still holds, and is 



a relationship between matrix elements in two different 
continuum theories, one unquenched and the other par- 
tially quenched. If one regularizes these two theories in 
the same way, then the relationship holds for all values of 
the renormalization scale [i. In particular, the anomalous 
dimension matrix of the four operators 09 ont should be 
the same as that of the original operators Oj . The results 
of Ref. [2] check this explicitly at one-loop order]^] 

It will be useful in the following to consider also matrix 
elements of 0? ont between kaons having different flavors 
and tastes. First we note that, because taste is a good 
symmetry in the continuum limit of the lattice theory, 
the matrix elements between type 1 and 2 kaons vanish 
unless they have taste P: 



(K° B1 \Of ont '\K° B2 )=0 if B^P 



(20) 



Here B labels one of the 16 choices of taste for the exter- 
nal kaons, as will be described shortly. Second, we con- 
sider matrix elements between mixed flavor kaons. Let 
K B12 be the kaon created by -Di^®^)^, and K B21 be 
the antikaon destroyed by D 2 (-j 5 <^^ B )Si. Here we are la- 
beling tastes by a hypercube vector B = (Bi,B 2 , B 3 , B4), 
in which each entry is either or 1, and 

e s = ef i ef 2 ef 3 €f% &. = tS). (21) 

Thus, for example, B — (1,1,1,1) corresponds to taste 
P. We then find that, for each value of j, 

^(K° P1 \Of ont '\K° P2 ), (22) 



(K° B21 \Of^'\K B12 ) 



where the sign s B is 



This result is obtained by Fierz transforming the 
operators in order to bring the bilinears into an 
"(SiD 2 )(S 2 Di)" form. One must simultaneously Fierz- 
transform in color, spin and taste. While the operators in 
Eqs. (9p2| are, by construction, Fierz- invariant in color 
and spin, they are not Fierz-invariant in taste. Taste P 
Fierz-transforms into all tastes, with the weight factor 
being s B /A. We stress again that the result (22) holds 
only in the continuum limit, for it relies on having an 
exact taste symmetry. 



C. Lattice operators for a/0 

Numerical calculations of the matrix elements required 
for Eq. ( 19 ) are being carried out in a lattice theory with 



three rooted sea quarks and two flavors each of valence 
down and strange quarks. This theory provides a lattice 
regularization of the PQ continuum theory described in 
the previous subsection. In this subsection we discuss the 
impact of the discretization errors inherent in the lattice 
regularization on the extraction of the desired matrix el- 
ements. The dominant issue is the presence of taste- 
symmetry breaking for a^O. 

To start with, we must choose a discretization of the 
continuum operators. The simplest choice is to use op- 
erators living on a 2 4 hypercube, using the method of 
Ref. to obtain operators with the desired spins and 
taste. We call the resulting operators 0^ at . The details 
of our particular implementation have been described in 
Ref. [14] and will not be important. What matters here 
is the structure of the matching between lattice opera- 
tors and those defined in the PQ continuum theory. The 
general form for the four-fermion operators i^] 



O 



ContA 



■o 



Lat 



+ — [taste P ops. 

47T 



4tt 



[other taste ops. 



-a [various taste ops.] + a [various taste ops. 



+ ..., (24) 

where the ellipsis indicates terms of higher order in ot and 
a. The symbol = means here the equality of the matrix 
elements of the operators on both sides of this equation, 
evaluated in their respective theories (PQ continuum on 
the left-hand side, lattice theory on the right). Thus all 
operators on the right-hand side are lattice four-fermion 
operators, and "taste B" indicates that both bilinears in 
the operator have this taste. The expression "various 
taste ops." implies that there are operators both with 
taste P and with other tastes. The set of operators which 
can appear is determined by the lattice symmetry group. 
We have separated out the one-loop contributions in 



sb = ^tr (£b£ 5 £b£5) • 



(23) 



(24) because they have been calculated in Ref. |14) . 



Strictly speaking, one must use a regularization which preserves 
the Fierz identities, such as the RI-MOM scheme of Refs. [231124] 
or the MS scheme proposed in Ref. 1251 . 



4 There are no contributions proportional to a because these would 
arise from mixing with dimension 7 operators, but, as explained 
in Ref. |15] . such operators have the wrong tastes to contribute 
to the desired matrix elements. 
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matching to a Fierz-invariant MS scheme in the contin- 
uum PQ theory. We use a/(An) (rather than just a) 
because the largest one-loop coefficients are of 0(1) x 
a/(47r). For the two-loop terms we do not know a pri- 



whether the corrections are 



*/(4^)] 5 



so 

we make the conservative choice and assume the former. 
Given the numerical value of a(l/a) for present lattice 
spacings (in, say, the MS scheme) , it is argued in Ref. [T3] 
that an appropriate phenomenological power counting is 
a /(An) ~ a 2 ~ (oAqcd) 2 <C 1. We adopt this power 
counting here, so all the displayed correction terms in 
Eq. p4| are formally of the same (small) size. 

We have further separated in ( |24| the one-loop con- 
tributions from operators with taste P from those from 
operators with other tastes. This is because, in our com- 
panion numerical calculations, we explicitly include (for 
practical reasons) only the taste P one-loop contribu- 
tions. In other words, the actual lattice four-fermion op- 
erator wc use is 



^-jLat, Actual ^jLat 



An 



[taste P ops.]. 



(25) 



Moving the other contributions in Eq. ( 24 1 from the lat 



tice to the continuum side of the equation (which can be 
done using tree-level matching as these contributions are 
of NLO due to the explicit factors of a/An, a 2 and a 2 ) 
we end up with 



.nLat, Actual ^^iCont' 



An 



[other taste ops.] 



-a 2 [various taste ops.] 
-a 2 [various taste ops.] ■ 



(26) 



Here operators to the right of the = are now continuum 
four-fermion operators. We see that our lattice operator 
corresponds in the PQ continuum theory to the operator 
we want together with several undesired operators. 
A similar analysis can be done for the bilinear opera- 



tors appearing in the denominator of Eq. ( 19 ). This case 



is simpler because, to all orders in perturbation theory, 
there is no mixing with other bilinears, due to the lattice 
symmetries |27| . Again, in practice we use a one-loop 
corrected operator, which can be written (for k = 1, 2) 



must use lattice symmetries, and appropriately general- 
ize the analysis given for the Bk operator in Ref. |15j . 
This exercise turns out, however, to be unnecessary when 
considering the Bj at NLO in SChPT. To explain this 
conclusion we must turn to the issue of mapping opera- 
tors into the chiral effective theory. 



III. MAPPING OPERATORS INTO SCHPT 

A. Review of SChPT 

We begin with a brief review of the relevant aspects of 
SChPT. More details are given in Refs. [13] and [30]. It is 
an effective theory constructed in three steps. First, one 
determines the Symanzik continuum effective Lagrangian 
describing the interactions of quarks and gluons with 
p <C I/a, which incorporates the leading discretization 
errors proportional to a 2 . Second, one maps the resulting 
theory into its chiral counterpart, in which the degrees of 
freedom are the pseudo-Goldstone particles produced by 
spontaneous chiral symmetry breaking. It is straightfor- 
ward to do this mapping only for an unrooted theory, i.e. 
one in which one keeps all tastes as dynamical degrees of 
freedom. The final stage is to account for the rooting of 
the quark determinant by including appropriate factors 
of 1/4 by hand for diagrams containing sea-quark loops. 
This last stage has been put on a firm theoretical footing 
by the work of Refs. [SI [32]. 



The standard power counting in SChPT is p l ~ m ~ 
a 2 . Here a 2 = a 2 av(n/a) 2 is the size of the leading taste- 
breaking corrections with HYP, asqtad or HISQ fcrmions. 
As described above, when one considers matrix elements 
one must also include taste-conserving discretization er- 
rors proportional to a 2 (without factors of a since HYP 
fermions are not fully improved) and matching errors 
proportional to a /(An) and a 2 . In the extended power 
counting introduced in Ref. [15] one assumes 



P 



a 
An 



I Lat, Actual 



[Sft(75®£5)A 



I Cont' 



(28) 



and 



We stress that the peculiar-looking choices a 2 <- 
a /(An) ~ a 2 are particular to the case at hand and are 



2 [n / ^, c \ n 1 2r „ • _ „ _ i /nn\ phenomenologically based. The choice a 2 ~ a 2 is made 

-a c [b k (75 ® &)D k \-a various taste ops. (27) £ ° ■ f ? 



with c an unknown constant of 0(1). In this case there 
are no errors proportional to a /(An). 

It is straightforward, although tedious, to enumerate 
the operators which appear in Eqs. ( [26] ) and (27) in the 
terms proportional to a /(An), a 2 and cP. For the a /(An) 
terms, the full list has been given in Ref. [25], along with 
their one-loop coefficients^] For the other operators, one 



5 The coefficients are given in Ref. 1281 only for the Wilson gauge 
action, rather than for the improved Symanzik gauge action used 
in practice. The results differ little, however [29] . In particular, 
the same operators have the largest coefficients in both cases. 



because it is found that taste-breaking discretization er- 
rors are numerically enhanced, and only after suppression 
by a 2 are they com parab le to other discretization errors. 
As explained in Sec. II C the choice a/ (An) ~ a 2 is based 
on the explicit results for one-loop matching coefficients. 

The Symanzik continuum theory obtained in the first 
of the steps described above is a partially quenched the- 
ory containing 28 quarks (3 sea and 4 valence, each with 4 
tastes) and 16 ghost quarks. It is convenient to collect the 
corresponding fields into a column- vector Q. In the com- 
bined chiral and continuum limit, the Symanzik action 
has a graded chiral symmetry, bU(28\16) L x bU(28\16) R . 
To display this we define left and right-handed Eu- 
clidean fields as usual, e.g. Ql = (1— 7s)/2Q and 
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Qr = Q(l-75)/2, so that 



^Sym 
The symmetry is 



QrPQr + QlPQl- 



(29) 



Q L -> LQ Ll Q R RQ R , Q l ->■ Q L tf, Q r ->■ Q R R\ 

(30) 

with L 7 R G 5?7(28|16)i ! fl. This graded symmetry is 
spontaneously broken down to its diagonal subgroup, 
leading to 44 2 — 1 pseudo-Goldstone particles^ 

The chiral effective theory contains only the light 
Goldstone particles that result after symmetry break- 
ing. These are collected as usual into a C/(28|16) matrix 
£ = exp(i$//) (with / such that f n w 132 MeV), where 



/ J7 7T+ if+ 

7T" D K° 

K- K° S 

\ : : : 



(31) 



Here each entry in the matrix is a 4 x 4 block correspond- 
ing to the 16 different tastes. Under the chiral symmetry, 
£ transforms as 



E — > LY.R). 
The LO chiral Lagrangian is 



(32) 



f 2 R f 2 
C x = y str (0 p E0 M £ t ) ^-str (ME + ME») 



+ ^[str(<I>)] 2 +a 2 (U+W), 



(33) 



where "str" stands for supertrace or "strace" , and M, is 
the mass matrix 

M = diag(m u ,m u ,m u ,m Ul m d ,m dl m d ,m dl m s , ■■■). 

(34) 

The mo term represents the effect of the axial anomaly 
(with normalization as in Ref. [13]); mo is to be sent to 
infinity to remove the unwanted non- Goldstone singlet 
particle [34] . The last term is the taste- breaking potential 
arising from discretization errors^ 

The taste-breaking potential consists of a single strace 



6 The grading does lead to some subtleties in the analysis of sym- 
metries and their implications, but these do not effect perturba- 
tive calculations in the resulting chiral theory 1331 . 

7 Although this contribution is proportional to a^, it is conven- 
tional to pull out an overall factor of just a 2 . This is a purely 
notational matter, since the difference can be absorbed in the 
LECs contained in U and U' . 



component IA and a double strace part W . The former is 
-U = C l str(^ n) E^Et N 



C 4 



2 E str (^ S ^ s+h - c - 

V 

+C 6 SstrtoE^Et), (35) 



where (following the notation of Ref. [TS]) £g is a 4nx4n 
matrix with the 4x4 taste matrix £b repeated along 
the diagonal n times. This potential contributes, along 
with the mass term, to pseudo-Goldstone particle masses, 
whose LO form (for flavor off-diagonal states) is exem- 
plified by 



L xy,B 



B (m x + my) + a 2 A(£ B )- 



(36) 



Here m xy> B is the mass of the pseudo-Goldstone bo- 
son composed of a valence quark of mass m x and a va- 
lence antiquark of mass m v , and having taste B. The 
taste-dependent discretization errors A(£b) depend on 
the LECs C\, C 3 , C4 and C 6 — explicit forms are given in 
Ref. [13] . In addition hi leads to four-pion vertices which 
contribute to the desired B— parameters at one-loop or- 
der. 

Both the mo term and the two strace potential W lead 
to "hairpin" (quark-line disconnected) vertices. Only the 
former contributes to the diagrams that enter here, and 
thus we do not reproduce the form of W . 

The procedure for accounting for rooting has been ex- 
plained in Ref. [13] . In essence, one must include by hand 
a factor of 1 /4 for each contribution which corresponds to 
a sea-quark loop. In the present calculation, it turns out 
that there are no diagrams containing sea-quark loops, 
as explained in Sec. IV A The only place where the 1/4's 
enter is in the quark loops implicitly contained in quark- 
line disconnected meson propagators. The impact of the 
1/4's is worked out in Ref. [15] . We quote the result 
only for the taste-singlet channel, since this is the only 
disconnected propagator we need. For a valence meson 
composed of a quark and its antiquark having mass m x 
converting to a similar meson composed of quark and 
antiquark of mass m y , the disconnected part of the prop- 
agator is (after sending mo — > 00) 



4 ( q 2 + U I )( q 2 + D I )(q 2 + S I ) 
'3 (q 2 + X!)(q 2 + Y!)(q 2 + ^){q 2 + m ) 



{37) 



Here we are using the same compact notation used to 
present results for B K in Ref. [TH]. Xj is the squared 
mass of the flavor off-diagonal, taste singlet pion created 
by .01(75 (g) ^/)Z?2, which at LO is 



Xi = 2B m x + a 2 A (6). 



(38) 



Yj, Uj, Dj and Sj are defined similarly. By contrast, 7r" 
and r\i are the mass eigenstates in the sea-quark sector, 
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which thus include hairpin contributions, and are given 
by 



Ut+Dj+Sj 1 



T ^Uf+Dj + Sj-UjDj-UjSj-DjS!, 

(39) 

with the upper sign for the ir® and the lower for the ijj. 
In the isospin limit m u = rrid one recovers the familiar 
results 7T° _ Uj and m _ ^ + 2S z )/3. 

The final issue to be discussed is the impact of using 
a mixed action, with different types of staggered valence 
and sea quarks. Here we can rely on the corresponding 
analysis for Bk |19j ■ The conclusion is that mixed-action 
effects can enter either through loops involving mixed 
valence-sea pions or through the presence of new hair- 
pin vertices of vector and axial taste. As will be seen 
in Sec. |IV A| however, in the present calculation there 
are, at one-loop, no contributions from mixed pions, and 
no contributions from vector and axial hairpins. Thus 
the only impact of using a mixed action is that the val- 
ues of the LECs associated with discretization errors are 
changed. This is not a concern, however, since these val- 
ues are to be determined by fits to simulation results. 



B. Operator mapping at leading order 



In this subsection we map the BSM operator s u sed in 
lattice calculations, i.e. the (^at.Actuai Q f Eq (|26|, into 
the chiral effective theory at leading order (LO). We must 
first map these operators into the Symanzik effective ac- 
tion. This is simplified by working at LO, which implies 
that we can drop corrections proportional to a/(4w), a 2 
and a 2 . It then follows from Eq. (26 1 that the LO map- 



ping is simply into the Of ont . Thus our task is to map 
the latter operators in the Symanzik effective theory into 
the chiral Lagrangian. 

The method for doing so was developed in Ref. pTj 
and used for the Bk operator in Ref. [T3]. One intro- 
duces spurion fields in such a way that quark-level op- 
erators become invariant under chiral transformations, 
then determines the LO operators in the chiral effective 
theory containing these spurions. Since all that matters 
are chiral transformation properties, the choice of color 
contraction is irrelevant, so both 2 ont and 0§° nt map 
into the same set of LO operators (of course with dif- 
ferent LECs). The same statement holds for 02° nt and 

£)Cont' 

It turns out to be simplest to map the component parts 
of the £>Cont' se p arate i y [ see Eqs. (9jl2l)]. We begin with 



O 



S+P 



2 ^Cont 



^Cont 



(40) 



where the factor of 2 is for later convenience. We do not 
specify the color contraction since the subsequent results 
hold for both choices. This operator is a component of 



2 ont and 0§ 0Dt . We first rewrite it in generic form 

O s +p = Qr{1 ® F 1L )Q L Q R (1 ® F 2L )Q L 

+Ql{1®F ir )Q r Q l {1®F 2R )Q r (41) 

where the second term on the right-hand side is the parity 
conjugate (p.c.) of the first. In the matrices (1 ® F) the 
"1" indicates the identity matrix in spin space, while the 
spurion F is a matrix in the tensor product of flavor and 
taste space. In order to obtain the desired operator, one 
must choose the spurions as follows: 

Fil,Fi R S^Sdij £5 and F 2 l,F 2R -> 5 uS2 8 D2t j £ 5 , 

(42) 

where i and j are flavor indices. 

Next we note that Os+p is invariant under chiral 
transformations if the spurions transform as 



F kL ^RF kL tf and F kR -> LF kR R\ 



(43) 



where k = 1,2. We now imagine that Os+p is inserted 
into the Symanzik action. Then the desired matrix ele- 
ments can be obtained by taking (a sum of) functional 
derivatives with respect to appropriate elements of the 
F k L and F kR and then setting the spurions to zero. 

The final step is to add to the chiral Lagrangian all 
operators composed of the new spurions together with 
S, £t, M, M\ derivatives and spurions coming from 
the 0(a 2 ) terms in the Symanzik action, such that the 
overall operator is invariant under chiral transformations, 
Euclidean rotations, and parity. Furthermore, we need 
only keep operators proportional to Fh j F 2 l and F\ R F 2R , 
since only these will survive when we take the functional 
derivatives of the chiral partition function needed to ob- 
tain the desired matrix elements. We see immediately 
that the LO operators will be those involving E and £T 
alone, with no derivatives, mass terms or a 2 spurions. 
There are two such operators: 



Ol = str(EF 1L £F 2L )+p.c. 



(44) 



where the parity conjugate is obtained by L O R and 
£ O £t and 



Ol = 



str (EF 1L )str(E J F 1 2L )+p.c. 



(45) 



One now sets the spurions to their original values, as in 
Eq. (42). The resulting two operators will appear with 



independent, unknown coefficients. 



To map 2 ont ' and O^ ont ' into the chiral theory we 
also need to consider Oj? nt . It turns out that this oper- 
ator maps into the same two chiral operators as Os+p- 
To see this, we note that, in addition to the form (151, 
the operator can be written 



TCont' 



/i7f75 ® #2(7^75 ® ^5)^2- 

(46) 

Combining these two forms one finds that the spurion 
representation of the operator is 



O 



T° nt ' = Qpilnlv ® F 1L )Q L Q R (w v ® F 2L )Q L + p.c. 

(47) 
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From the point of view of chiral symmetry, this operator 
transforms in exactly the same way as Os+p, Eq. (41 1. 



Thus its mapping into chiral operators has the same form. 

The final stage of the mapping is to note that the rela- 
tive coefficient of the two chiral operators 0* and O x is 
fixed, so that there is only one overall LEC. This holds 
only for the particular linear combinations of Os+p and 
Oj. ont ' that appear in 0§^ nt ' ■ The key observation is that 
the coefficients of the two chiral operators are exactly the 
same as they would be if one set a = in the Symanzik 
theory. This is because all factors of a are explicit, and 
there are none in either chiral operator. But setting a = 
in the Symanzik theory leads to the PQ continuum the- 
ory considered in Sec.|II B Since taste symmetry is exact 



in this theory, the results ( 20 ) and ( 22 ) must hold for the 



matrix elements at both the quark level and the chiral 
level. Furthermore, they must hold order by order in 
the momentum-quark-mass expansion of SChPT, and in 
particular they must hold at LO in the standard power 
counting. The leading-order matrix elements are simple 
to calculate. If we write the chiral mapping of O^ 011 * as 



0* = (c 2a O x + c 2b O x ) , 



(48) 



then we find 



(K m \O x \K B2 ) ho 
(T? B21 \O%\K B12 } LO 



32 

F 

8 
J 



-C2a-M£,B^B^)- (50) 



These results are consistent with (20 1 and (22) only if 
c 2 f,. We thus conclude that 



C-2u 



(51) 



The same form holds for Og but with a different coeffi- 
cient, C3. 

The fact that, at LO, there is only one unknown LEC 
could have been anticipated from the result that there is 
only a single LEC in the mapping of O2 into continuum 
ChPT [16] . We also note that a similar analysis holds for 
the chiral mapping of the Bk operator |15j . 

We now turn to the chiral mapping of O4 | nt . These 
are composed of 



S -p = 2 O 



^Cont' 



-o 



Cont' 



and 



O 



V-A 



= 2 



0£ ont ' - o A 



Cont' 



(52) 



(53) 



In terms of spurions, the former operator is 

Os-p = Qr{1 ® F 1L )Q L Q L (1 ® F 2 r)Qr + p.c. (54) 

Here, the spurions F±l, F2R and their parity conjugates 
transform as above [E q. ([43] )] , and are set at the end to 
the same values as in (1421). Note, however, that Eq. (54) 



differs from the spurion form of Os+p, Eq. (41). The for- 



mer is proportional to F^Fr, while the latter to F^F^. 
This leads to the presence of only a single LO chiral op- 
erator, 

O^ = str (EFi L ) str (p F 2 r) + p.c. (55) 
Turning now to the V — A operator, its spurion form 



is 

Ov-A 



J2Ql{1v® F 1L )Q L Q fl (7 M ® F 2R )Qr + p.c, 

(56) 

where the spurions now transform as 

F kL ^LF kL tf and F kR -> RF kR rf . (57) 

At the end they are set to the same values as the other 
spurions, 

Fil,Fir 5i t siSm,j £5 and F 2 l,F 2 r -> £,„s2<5D2,j £5- 

(58) 

The single LO chiral operator that this maps to is 

O x d = str (tfF 1L ZF 2R ) + p.c. (59) 
Combining these two operators into O4 g nt and enforc- 



(49) ing the relations (20) and (22) from taste symmetry and 



Fierzing, we find again that the coefficients are related. 
The chiral mapping is to 



^4,5 



C4,5 (Oi 



o x d ) 



(60) 



Finally, we need to map the pseudoscalar densities ap- 
pearing in the denominator of Eq. ( 19 1 into the chiral 



theory. Since we are working at LO, the operators in 
the Symanzik effective theory that we need to map are 
[Sk(l5 ® ^5)F> k } Cont ' for k = 1,2. This is a standard 
exercise and we find 



Cont' 



U: 



(61) 



0* u = str(SF feL )-p.c. (62) 

Expressing the constant c^q in terms of other LECs is 
not useful here since the corresponding constants in the 
mapping of the numerator of ( 19 1 are unknown. 



C. Mapping at next-to- leading order 



An important conclusion from the previous subsection 
is that, because the LO chiral operators Oj contain no 
derivatives, they give rise to LO matrix elements that are 
non- vanishing in the chiral limit. This is seen explicitly in 
Eqs. ((49]) and (50). Unlike for the operator, there is 



no chiral suppression, a result that is well known in con- 
tinuum phenomenology. This means that higher-order 
chiral operators which contain factors of M, a 2 , a /(in), 
a 2 , or which contain derivatives can only give rise to ana- 
lytic contributions to the Bj . Non-analytic contributions 
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at NLO can arise only from one-loop diagrams involving 
the LO chiral operators. 



IV. SCHPT RESULTS FOR B PARAMETERS 
AT NLO 



Because of this, we can determine the functional form 
of the NLO contributions to the Bj without explicitly 
enumerating all the higher-order chiral operators which 
appear when we map the 0^ at ' Actual into the chiral the- 
ory. We know that operators which come with two 
derivatives will lead to analytic terms oc m? K , while op- 
erators arising from discretization errors in the action or 
the operators lead to analytic terms oc a 2 . This holds for 
both the numerators and denominators of the Bj. 

We can also see that there are no NLO analytic terms 
oc a/(47r). These only arise from the "other taste" oper- 



ators of Eq. (26), which enter in the numerators of the 



Bj. When we match these four-fermion operators from 
the lattice to the Symanzik effective theory, the LO con- 
tinuum operators that result will still have tastes other 
than P (since taste-breaking would bring in further fac- 
tors of a 2 ). Thus their matrix elements with taste P ex- 
ternal kaons vanish at LO. This in turn implies that the 
LO chiral representation of these operators (which con- 
tain no derivatives) must have vanishing tree-level ma- 
trix elements between kaons of taste P. Their one-loop 
matrix elements will be non- vanishing, but, because of 
the overall factor of a/(4ir), these contributions enter at 
next-to-next-to-leading order. 

Finally, we consider terms involving insertions of the 
quark mass matrix M. Here an explicit enumeration is 
useful. We find two types of chiral operators. First, those 
in which the LO operators are multiplied by str(M£i) 4- 
p.c. These lead to analytic corrections oc m u + + m s . 
Second, factors of M£t an d EM* can be inserted in the 
LO chiral operators. The corrections to 0^ 3 are, for 
example, 



str (Ei^EFsiEM 1 ) + (1 O 2) + p.c. (63) 
str (EFxiEMt) str (£F 2L ) + (1 o 2) + p.c. (64) 



The first of these operators does not contribute to the 
desired matrix element at tree-level, while the second 
gives a contribution oc m x + m y . There arc no terms 

oc m x 



' y 1 



as can also be seen more generally using 
CPS symmetry [35] . A similar analysis leads to the same 
conclusion for the corrections to 0\ 5 . 

In summary, NLO analytic corrections to both the nu- 
merators and denominators of the Bj are proportional 
to m 2 K , a 2 , a 2 a , a 2 , m u + m d + m s and m x + m y . 
Since we use pseudo-Goldstone external kaons, for which 
m 2 K oc m x + m y at LO (with no a 2 terms), the m 2 K and 
m x + m y corrections can be combined into a single term. 



The analysis of the previous section shows that, to 
NLO in SChPT, we have 



B< = 



(K P1 \Of\K P2 )^ loop 



3 Nj (KUo^lOh-iooAOn^l 
+ analytic NLO, 



L P2/ 1— loop 



(65) 

with the chiral operators in the numerator defined in 
Eqs. (51 1 and (60), while those in the denominator are 
given in Eq. (62). The subscript "1-loop" here means 



the sum of tree-level and one-loop contributions. At LO, 
the matrix elements in both numerator and denomina- 
tor are constants, independent of quark masses and kaon 
momenta. Explicitly, we find 



LO= 2 Tc.32// 2 
Nj (8*c bil //) 2 



= ± 



Njc 2 n 



(66) 



Throughout this section the upper sign applies for j = 
2, 3, and the lower sign, for j = 4, 5. The result (66) has 



no predictive power since we do not know the constants 
Cj. The B}j° are simply the values of the Bj in the joint 
chiral-continuum limit. 



The predictive power of Eq. (65) arises because the 



one-loop contributions involve only the LO chiral opera- 
tors, implying that the relative contribution of the chiral 
logarithms is determined. 

The tree-level and one-loop contributions to the kaon 
matrix elements in the numerator of Bj are shown in 
Fig. [T] Here we distinguish between the single and double 
strace sub-operators contained in the Oj. For example, 



from Eq. (51 ) we see that the single strace component of 



0£ is Ca °f Eq. ( 44 ) while the double strace component 



is Oi of Eq. (45 1. Given the tastes of the external kaons, 



only the double strace components contrib ute at tree- 
level, as shown for the case of 2 by Eq. (50). Both 
components contribute at one-loop order, as shown in 

The contributions to the denominator are shown in 
Fig. [2j They are simpler since there is only a single strace 
component. Noting that the square boxes in Figs. [T] and [2] 
correspond to identical chiral operators, and accounting 
for the fact that the loops in Figs.JTJb) and (c) can appear 
on either external kaon line, we see that the contributions 
to Bj from Figs, [ljb) and (c) cancel with those from 
Figs.j2jb) and (c). Wavefunction renormalization factors 
also cancel. These are the same cancellations as those 
found for Bk in Ref. [TS]. Thus we need only keep the 
diagrams of Figs. [TJd) , (e) and (f). 

It is useful to draw the quark-line diagrams contribut- 
ing to Figs, [ljd-f). These are shown in Fig. [3} We recall 
that these are primarily a device for tracking the flavor 
indices of mesons in the diagrams that contribute to the 
SChPT calculation. They also correspond, however, to 
different ways of routing the quark propagators of the un- 
derlying lattice calculation so as to make loop diagrams. 



10 



(a) 



(b) 




-Q- 



O 

(f) 



(e) 



FIG. 1. SChPT diagrams contributing to the numerator of 
Bk ■ (a) tree-level; (b-f) one-loop. The double strace compo- 
nents of the chiral operators are represented by two square 
boxes, one per strace, while the single strace components are 
shown with one rectangular box. (It turns out that there is 
no contribution from a diagram of the form of (f ) but with 
a two strace operator.) The filled circle is the full LO vertex 
from the SChPT Lagrangian, including 0(a 2 ) terms. For (b) 
and (c) we have not shown separately the diagrams in which 
the loop is attached to the other external kaon. 



(a) 



-QQm 

(b) (c) 



FIG. 2. SChPT diagrams contributing to the single kaon 
matrix elements in the denominator of Bk'- (a) tree- level; 
(b-c) one-loop. Notation is as in Fig. [I] 




FIG. 3. Quark-line diagrams contributing to (a) Fig. [TJd), 
(b) Fig. [lje), and (c) Fig. |TJf ) . In each case, each of the 
filled boxes correspond to one of the factors of E (or E*) in 
the chiral operators. The square boxes arise from the double 
strace operators and thus have flavor SiDi or 52-D2- The 
octagonal boxes arise from the single strace operators and 
thus have flavor S1D2 or S2-D1. 



having different flavors (S2D2 versus SiDx). Further- 
more, because both the external kaon and the bilinear 
represented by the square box have taste P, the meson 
in the loop must be a taste singlet. An important corol- 
lary is that the second complication due to the use of a 
mixed-action — namely, the change in taste- V and taste- 
A hairpin vertices — does not impact the present calcula- 
tion. 



In the latter interpretation, each of the boxes corresponds 
to one of the component bilinears in the four-fermion 
operator, and in the case where the boxes are octagons 
rather than squares one must first Fierz transform the 
operator into its (SiZ^XS^-Di) form. 

We see from Fig. [3] that there are no diagrams in- 
volving valence-sea mesons. Such mesons do contribute 
to Figs, [ijb) and (c) and to Figs. [2^b) and (c), but 
these contributions cancel at NLO as discussed above. 
This is an important simplification because it means that 
we do not need to determine the masses of valence-sea 
mesons, which contain a different a 2 contribution from 
the valence-valence mesons. Thus one of the potential 
complications from using a mixed action does not occur. 

We also see, from Fig. |3ja), that Fig. [ijd) involves 
only the quark-disconnected, hairpin part of the meson 
propagator. This holds because one cannot have a quark- 
connected propagator joining mesons composed of quarks 



A. NLO SU(3) SChPT result 

In this section we give the general form of the next-to- 
lcading order corrections for a 1 + 1 + 1 flavor theory, i.e. 
the rooted theory in which we keep the sea quark masses 
general. 

We break up the corrections as follows 



B 3 = Bj° [l + SB 



anal 



6B c 3 onn 



SBf lsc ] 



(67) 



where the first correction contains the analytic term, 
the second is the contribution from Figs, [lje) and (f), 
which we refer to as "connected" since they do not in- 
volve hairpin vertices, and the third is the contribution 
from Fig.JTJd), which is "disconnected" as it involves only 
hairpin vertices. 



The analytic terms have been discussed in Sec. Ill C 
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and have the form 

SBj nal = Cji(m x + my) + c j2 (m u + m d + m s ) 

+c j3 a 2 + c j4: a 2 a + c j5 a 2 . 

We find the connected contributions to be 
-1 1 



(68) 



SB C j° nn 



(4tt/)2 16 ^ 
±£{X B )±£{Y B ) 



[2£(K B ) - 2K P £(K B ) 



(69) 



where B is a taste label, which is summed over the 16 
possibilities. We use the following abbreviations for me- 
son mass squareds [19]: X B = Tn%. x B , Y B = m y B and 
K B = m 2 y B . The chiral logarithmic functions are, in 
infinite volume, 



£{X) =Xln 



X 

J 2 



£(X) 



In 



X 

7 2 



i. 



(70) 



where [i is the renormalization scale in dimensional reg- 
ularization. The dependence on fi is absorbed by the 
implicit /i dependence of Cji and (for j = 4,5) Cj±. Fi- 
nite volume corrections to these logarithms are standard 
and are given, e.g., in Ref. |15j . 

We find the disconnected contributions to be 



SBf sc = ±~ 

3 4 J2 



W) 4 



[DL(q) + D I yy (q)+2Dl v (q)}, 



where the disconnected propagators are given in Eq. 
ft is straightforward, though tedious, to evaluate these 
integrals by summing the contributions from the poles in 
the propagators, following the method of Ref. [33]. For 
the sake of brevity, however, we only quote the result in 
the isospin symmetric limit (i.e. for a 2 + 1 flavor theory): 

' ;{I V +Ix 



SB ? sc=± 3(^fr 



Iy) 



(72) 



where 



I v = £{ m ){L I - m ){S I - 



■Vi) 



X, 



Y, 



~Vi 



(73) 



and 



ix 



{L,-X I ){S I -X I ) 
(Vi-Xj) 

2 1 



£{Xj) + £{Xi)x (74) 



1 



1 



Vi 



-x, 



(75) 



Yj—Xj Lj — Xj Sj — Xj 

I Y =Ix{X^Y). 

Despite appearances, this expression is finite when Xj — > 
r/j, Lj, Sj or Yj, and similarly for Yj. The ln/i depen- 
dence is proportional to 

Li + Si-th-Xi- Yi 
= h$L t + S I )-{Xi+Yi) 



(76) 



B 



(m u +m d +m s ) - (m x +m y ) 



-a 2 A(6), (77) 



and can thus be absorbed by shifts in Cj%, Cji and 0,-4. 

We can check our result with that of Ref. [TB] by taking 
the continuum, unquenched limit. In this limit, all tastes 
are degenerate, and 



DIM 
Dly(q) 



-2 



1 



1 



q 2 + m\ 



3<7 2 



1 



3 q 2 + m 2 
1 



q 2 + S 3q 2 +m 2 ' 



(78) 
(79) 
(80) 



Thus 



5Bf nn 



while 



(4tt/)2 



2£(K) - 2K£(K) ± £{m\ ) ± £{S) 



5Bf isc 



■(47r/)s 



.^( m 2)-f(S) + ^( m 2) 



(81) 



(82) 



Combining these results one finds that the unphysical 
logarithms £(S) cancel, and the results agree with those 
of Ref. [IS]. 



B. SU(2) SChPT 

The utility of SU(3) ChPT at the physical kaon mass is 
unclear, given the relatively large value of the expansion 
parameter m 2 K /K 2 i (A x <~ 1 GeV). Thus it has proved 
very useful to consider the strange quark as heavy and 
work instead with SU(2) ChPT. In this limit only pion 
loops lead to chiral logarithms, while the kaon is a static 
source. SU(2) ChPT was first developed in the contin- 
uum in Ref. [36], and extended and applied to Bk (and 
other quantities) in Ref. [37] . The extension to staggered 
ChPT was described in Ref. [12] in the context of the ap- 
plication to Bk- Our present application is very similar 
to (and indeed somewhat simpler than) Bk , and thus we 
can take over much of the work of Ref. [T5] . 

We first make clear the limit that we are considering 
in our extended partially quenched set-up. We take the 
valence and sea strange-quark masses (m v and m s , re- 
spectively) to be heavy, while m x , m u and remain 
light. All particles containing one or two strange quarks 
of either type are treated as heavy, i.e. the valence and 
sea quark kaons, yy particles and the rjj. The masses of 
these heavy particles are considered to be of the same 
order as A x . 

It is argued in Ref. [lJ5] that for Bk one can obtain the 
general NLO SU(2) SChPT result from the NLO SU(3) 
SChPT result using the following recipe: take the limit 
TOtt <C rriK and treat (m^/rriK) 2 as a small parameter 
of size (p/A x ) 2 , drop chiral logarithms of heavy parti- 
cles, and replace all LECs with (unknown) functions of 
m y and m s . The only loops that remain are those of 
mesons containing light quarks alone, e.g. xx particles 
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and sea-quark pionsj^] The argument holds also for the 
present application. Indeed, the argumentation is sim- 
pler here because the LO matrix elements are not chi- 
rally suppressed. In light of this simplicity, and because 
the treatment of Ref. [T5] suppressed some details, we 
explain the argument here. 

First, however, we give the result, as this will facilitate 
the subsequent explanation. Here, the full 1 + 1 + 1 flavor 
result is simple enough to present. It has the same form 
as the SU(3) result, Eq. @, 



LO,SU(2) 



1 + SB 



anal,SU(2) 



-SB 



conn,SU(2) 



disc,SU(2) 



SB 



but now the overall constant, £?lo,su(2)^ ^ g vame £ 

Bj in the SU(2) chiral limit (and for a = 0) rather than 

the SU(3) limit which applies for Bj°. ^ LO ' SU(2) thus 
has an unknown and general dependence on m y and m s . 
The analytic term becomes 



(83) 



SB 



anal,SU(2 



' = djirrix + d j2 (m u + m d ) 



+dj 3 a 



IJ4.C 



d j5 a 



(84) 



with new LECs. Because only the light meson loop sur- 
vives, the connected contribution simplifies to 



£gConn,SU(2) _ _|_ 



(i^W- ,85, 



Here f 2 is the decay constant in the SU(2) chiral limit. 
Finally, the disconnected contribution can be obtained 
from Eq. ( 71 ) by noting that the light-particle poles in 



D xy come with residues having an additional suppres- 
sion factor m Ut d/ m s,yi while those in Dyy (which are only 
present in a 1 + 1 + 1 flavor theory) are even more sup- 
pressed. Thus only the D\. x term contributes, and using 
S 1 [rf = 3/2 and 2ttj = XJj + Dj [which follow from 
Eq. (39) for m s ^> m u ( j\, we find 



SB 



disc,SU(2) 



± 



£(Xj) 



{U I -X I ){D I -X I ) 



2{^f 2 ) 2 \ y " 1J (ttz-Xx) 

(7T/-A/) 2 



(86) 



This expression simplifies in the isospin symmetric limit 
to 



SBf sc ' SU( - 2) (m u =m d ) 



± 



2(4tt/ 2 )2 



I (Xi) + (Dj - Xj)£(Xj) 



(87) 



There are two striking features of these results. First, 
the chiral logarithms for B 2 and B3 have exactly the 
same form but opposite sign to those for B4 and B$. 
This is not true for the SU(3) result. Second, the chiral 
logarithms for B 2 and B3 are identical to those found for 
Bk in Ref. [TH] Q These features will be explained by the 
following analysis. 

We now turn to the arguments which justify the pre- 
scription used to obtain the above results. The assump- 
tion made in Ref. Q23 is that the NLO SU(2) SChPT 
results can be obtained by working in SU(3) SChPT to 
all orders in m y and m s while working only to NLO in 
the small quantities ra x , m u , rrid, a 2 , a 2 and a 2 . This 
assumes that there are no non-perturbative contributions 
in m y and m s . In this approach, any number of loops of 
particles containing strange quarks are allowed, but only 
a single loop containing light particles. One then imag- 
ines taking the SU(2) limit, so that heavy particle loops 
lead only to contact terms, which are analytic in the light 
quark masses. 

It is important to note that, while loops of heavy parti- 
cles are not suppressed (since they give contributions like 
(m K /A x ) 2 ln(m^-//x) ~ 1), they are also not enhanced in 
the SU(2) power counting. Thus a contribution to the 
NLO SU(3) SChPT result whose suppression is by one of 
the small quantities in SU(2) SChPT remains at most of 
NLO in SU(2) power counting with the addition of any 
number of heavy loops. This implies that contributions 
of NLO in SU(2) power counting can arise either from a 
light particle loop or from an operator bringing in an ex- 



plicit factor of m x , m u , rrid, 



or a (in both cases, 



with any number of heavy loops) , but not from both. For 
the analytic contributions to Bj, this leads immediately 



to the result of Eq. (84), with dji — dj% being arbitrary 



analytic functions of m y and m s . 

A further observation is that, at NLO in SU(2) power 
counting, where heavy loops are collapsed to a point, all 
one-loop diagrams are tadpoles, i.e. with no vertices on 
the light particle propagator. They have the generic form 
of Figs, [ljd) or (e) with a pion in the loop. Diagrams of 
the type of Fig. [IJf) collapse to an analytic contribu- 
tion since there is necessarily a kaon in the loop. This 
observation implies that, since the loop itself gives a con- 
tribution of NLO in SU(2) power counting, the vertex 
to which it attaches must be LO in this power counting. 



t\b particles for B/J are light in the SU(2) limit, but these do 
not appear in the expressions for the Bj. The same is true for 
mixed mesons composed of a light sea quark and a light valence 
antiquark, or vice-versa. 



This holds for the full 1 + 1 + 1 theory, although the results of 
Ref. I19| enable a check only for the 2 + 1 flavor theory. 
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Since at LO there are no factors of a 2 , a 2 or a 2 in the op- 
erators, the vertex must be the same as that obtained by 
taking the SU(2) limit of the continuum result obtained 
to all orders in SU(3) ChPT. By assumption, however, 
this limit gives the same vertices as directly working with 
continuum SU(2) ChPT. Thus we conclude that the non- 
analytic NLO contribution can be obtained by calculat- 
ing the tadpole vertices using continuum SU(2) ChPT 
and inserting the propagators from the SU(2) limit of 
SChPT. 

The final step of the argument is to note that in SU(2) 
ChPT the relative contribution of the chiral logarithms is 
independent of the masses m y and m s . This point, de- 
rived in Ref. [37], will be explained briefly below. This 
relative contribution is therefore the same as in a theory 
with m x ~ m u ~ <C m y ~ m s <C Aqcd- In such a 
theory it is justified to obtain the NLO SU(2) SChPT re- 
sult by calculating at NLO using SU(3) SChPT and tak- 
ing the SU(2) limit, which is exactly the procedure used 
above except that the LECs are independent of m y and 
m s . We therefore conclude that the form of the SU(2) 
SChPT chiral logarithms obtained above is correct. We 
stress again that the LEC f 2 does depend on m y and m s ; 
what is independent of these masses is the remainder of 
the expressions in Eqs. ( 85p7 ). 



usual, we couple K into the chiral theory by introducing 
u = VS = cxp[i$/(2/2)], which transforms as 



In order to check the final step of this argument, we 
have calculated the chiral logarithms in the hybrid theory 
described above, in which the vertices are from contin- 
uum SU(2) ChPT while the propagators are from SChPT 
in the SU(2) limit. This calculation is also useful because 
it shows why the results for B 2 z (-84,5) have exactly the 
same (opposite) form to those for Bk ■ 

The set-up for continuum SU(2) ChPT in our case re- 
quires a generalization of the method of Ref. [37] , because 
our continuum theory is the enlarged, partially quenched 
one described in Sec. |IIB] In particular, we must account 
for the extra taste degree of freedom, the two types of va- 
lence S and D quarks, and the fact that our operators 
have a different form [Eqs. (9]fl2 1 compared to Eqs. (JlJ- 
[4J]. These features imply that the kaon fields can be 
collected into a rectangular matrix, K ab ~ £ a Sb, where 
the first index runs over the flavors and tastes of light 
quarks and ghosts, while the second runs over the two 
types of strange quark and their tastes. Thus a takes 24 
values (valence down quarks of both types, correspond- 
ing ghosts, and up and down sea quarks, all with four 
tastes) while b takes 8. The approximate chiral symme- 
try is now SU(16\8)l x SU(16\8) r . We only know how 
K transforms under the vector subgroup, and also under 
the SU(8) group rotating between the 8 different strange 
quarks: 

K — > UKV\ u e SU(16\8) V , V e 51/(8). (88) 

The construction of the SU(2) chiral Lagrangian for pions 
follows exactly the same steps as described in Sec. |III A[ 
aside from the absence of a 2 terms. Thus S is now a 



Lutf = UuB). 



(89) 



Thus the combinations uK and wK transform simply 
under the full chiral and SU(8) groups 

[uK) -> L{uK)V\ {u^K) -> R{u^K)V^. (90) 

The leading order kaon Lagrangian is then 

C X , K = tr {DptfDpK + m 2 K K^K) , (91) 

where the trace is over the SU(8) indices (which are not 
graded). The covariant derivative involves the u fields 
and is defined, e.g., in Ref. [37] , 

To map operators into this SU(2) chiral theory we use 
a variant of the spurion method described earlier. For 
example, we write one of the component operators of 
0§ ont ' as 

O s+ p = S(l <g> F 1L )Q L S(1 ® F 2L )Q L + p.c, (92) 



to be compared with Eq. (41 1. Note the absence of a 



chirality subscript on the heavy S field. For economy of 
notation, we reuse the symbols F\l, Ql etc., although 
they have different meanings in the SU(2) theory. Here 
S is a row vector containing the 8 strange quark fields, 
while Ql is a column vector containing the 24 light left- 
handed fields and F\l and F 2 l are 8 x 24 matrices. To 
obtain the desired operator the spurions have to be cho- 
sen as in Eq. ( 42 1 , where now the Kronecker 5 is a rect- 



angular matrix. Under chiral transformations, the light 
quark fields transform as in Eq. (30), while under SU(8) 
transformations S — > SV^ . Thus the spurions here must 
transform as 

F kL ^VF kL L? and F kR -> VF kR B? (93) 

for k = 1,2. The mapping into the chiral theory follows 
the same logic as for the SU(3) theory. The simplest 
operators that are allowed are 

x,su{2) = tT ( FlLuK )tr(F 2L uK) 

+tv(F 1R ^K)tr(F 2R u^K) (94) 



and 



X,SU(2) 



ti(F 1L uKF 2L uK) + ti(F 1R u^KF 2R u^K). 



O 

(95) 

By the same reasoning as before, 0^f ont also maps into 
the same two operators. Furthermore, calculating the 
LO matri x ele ments between kaon states and enforcing 
(Eol and Q, one finds that the full operator 0§ ont 



24 x 24 graded matrix, transforming as in Eq. ( 32 ) . As 



maps into the sum of the two chiral operators. There is 
thus only a single overall LEC. The same form holds for 

In SU(2) ChPT, operators with additional covariant 
derivatives acting on the kaon fields are not suppressed. 
As explained in Ref. [37] , however, by using the equations 
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of motion one can reduce these operators down to those 
without derivatives, up to contributions of higher order 
in the SU(2) power counting. This is the result that 
shows how arbitrary powers of m K can appear without 
impacting the coupling to pions. 

We now sketch how the mapping changes for 0^ nt . 
Here the quark-level operators, shown in Eqs. (54 1 and 



(56), contain both a Ql and a Qr field, in contrast to 
the LL or RR structure that appears for 0§ | nt ■ The 
net result is that one must build the chiral operators out 
of Fiz and F 2 r (or their parity conjugates). Note that, 
since Sl and Sr transform in the same way, we do not 
need to introduce new spurions F)~l and F^r- We then 
find that both Os-p and Oy-A are mapped into a linear 
combination of 



x,su(2) = tv ( FlLuK ) t ^ F2RV ) K ) + p. c . 



and 



O 



X,SU(2) 



tx(F 1L uKF2Ru'K) + p.c. 



Enforcing <|20h and (|22h we again find that Of ont ' 



4,5 



(96) 



(97) 



both 



map into the sum of these two chiral operators. 

The mapping of the bilinears in the denominator of the 



Bj is simpler and leads to [cf. Eq. (62)] 



tv(F kL uK) + p.c. 



(98) 



Now we observe that the form of the SU(2) ChPT op- 
erator for C^ nt ', i.e. the sum Of SV{2) +Of SU{2 \ is 
identical to the operator which represents the operator 
appearing in the numerator of Bk, and whose explicit 
form is given in Ref. [19]. This is because all three oper- 
ators have a LL+RR structure in terms of the light quark 
fields. Thus the one-loop corrections to the numerator 
of i?2,3 are identical to those for the numerator of Bk- 
As for the denominators, they have a somewhat differ- 
ent form (the denominator for Bk is mapped into SU(2) 
ChPT in Ref. [37]), but at NLO, and when expressed in 
terms of $ and K fields, they are proportional. This im- 
plies that the chiral logarithms in both numerator and 
denominator will be the same for B 2 , B% and Bk- This 
holds also in the hybrid calculation using SChPT prop- 
agators, since the vertices and the propagators are the 
same in all cases. We thus can take the SU(2) SChPT 
result from Ref. [T!5] and indeed find, as observed above, 



that it agrees with the sum of Eqs. (85) and (87). This 
checks our argument based on taking the SU(2) imit of 
the SU(3) SChPT expression. 

Finally, we can now understand why the chiral loga- 
rithms in B4 5 are, in the SU(2) limit, exactly opposite 
to those for B 2 ,3 and Bk- The first step is to notice that 
the LR + RL structure of the quark-level operators in 
i?4,5 maps into operators at the chiral level with one it 
and one u* . This is in contrast to the two it's or two 
it's for B 2 z- Next we note that there is a significant 
cancellation of chiral logarithms between the numerators 
and denominators of the Bj. In particular, as observed 



for Bk in Ref. [37], the denominator cancels the contri- 
butions from the numerator in which one of the u/w's is 
expanded out to C($ 2 ), while the other is unity. (This 
type of contribution is only non-vanishing for the two 
trace chiral operators in the numerator.) Thus the only 
contributions surviving the cancellation come from terms 
in which both u/iv's are expanded to linear order in $. 
Recalling that it = exp[i$/(2/)], we see that the expan- 
sion of the u x operators will lead to the opposite sign 
to that from u x it or iv x w . This implies that the chiral 
logarithms are opposite for B 2 3 and -84,5. 



C. Continuum PQ results 

As noted in the introduction, the continuum PQChPT 
result is not available in the literature, either for SU(3) or 
SU(2) ChPT. We can obtain these results by taking the 
continuum limit of our general formulae. In this limit, all 
taste breakings vanish, so we can make the substitutions 

K B -> K, Xr — » X, Y B -> Y, etc., (99) 

as well as setting a and a to zero. 
For the SU(3) case, we find 

^anai ^ Cjl (rn x + m y ) + c j2 (m u +m d + m s ) (100) 



2£(K)-2Kl(K)±£(X)±£(Y) 



(101) 



These results hold for the general 1 
We stress that in the result for SB^° n ' 



1 + 1 flavor theory. 
, K is the squared 
mass of the partially quenched kaon (composed of a quark 
with mass m x and an antiquark of mass m y ) as opposed 
to the physical kaon composed of sea quarks. 

For 8B^ lsc , the result in the 2 + 1 flavor theory is iden- 
tical to that in SChPT, given in Eqs. (72 - 75), except 
that the subscript I is dropped. We do not quote the 
1 + 1 + 1 flavor result explicitly, as it is lengthy, but it 
can be obtained straightforwardly from Eq. ( [7l] ) using the 
propagators of Eq. (37) with the subscript / dropped. 



The result for partially quenched SU(2) continuum 
ChPT in the 1 + 1 + 1 flavor theory is 



± 



LO,SU(2) 



1 + djxm x + d j2 (m u + m d ) 



£(X) 



(U-X){D-X) 



2(4^/ 2 ) 2 { n-X 
(U-X)(D-X) 



AX) 

-£(*) 



(n-X) 2 
(U-tt)(D-tt) 



(102) 
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This reduces in the isospin limit (m u = rrid = me) to 



LO,SU(2) 



1 + djim x + d^lmi 



± 



2(4tt/ : 



\— {e(X)(n-X)-£(X)}. (103) 



V. CONCLUSIONS 

We have presented the next-to-leading order results in 
staggered chiral perturbation theory for i?-parameters of 
the kaon mixing operators that generically arise in mod- 
els of new physics. We have done so for both SU(3) and 
SU(2) chiral perturbation theory. These results can be 
used to extrapolate lattice data obtained using staggered 
fermions to the physical light and strange quark masses. 
As a side product, we also provide partially quenched re- 
sults for both SU(3) and SU(2) ChPT in the continuum. 

We find that the results are much simpler in SU(3) 
SChPT than for Bk- Terms induced by discretization 
and matching errors in the lattice operators enter only 
through analytic terms rather than through chiral log- 
arithms. We also find that the use of a mixed action 
does not change the form of the NLO results. For SU(2) 
SChPT the results are of comparable simplicity to those 
for Bfc. Indeed, the chiral logarithms for B 2 and B 3 are 



identical to those for Bk, while those for B^ and B§ are 
opposite. In both SU(3) and SU(2) SChPT, if one works 
at fixed lattice spacing, the NLO expressions have the 
same number of unknown constants as those in the con- 
tinuum, as long as one first determines the masses of the 
valence pions and kaons of all tastes. 

It was pointed out in Ref. [16] that certain combina- 
tions of B-parameters have vanishing or small chiral log- 
arithms. The former combinations, dubbed "golden" in 
Ref. [16], remain golden in SChPT. The two examples 
built from B-parameters alone are the ratios B 2 /B 3 and 
B4/B5. The "silver" combinations are (one of B2.3) x 
(one of -84,5), for which the SU(3) chiral logarithms 
largely cancel. These turn out to be golden in SU(2) 
SChPT. It may be useful to use these combinations to 
improve the chiral extrapolations. 



VI. ACKNOWLEDGMENTS 

The research of W. Lee is supported by the Creative 
Research Initiatives program (3348-20090015) of the 
NRF grant funded by the Korean government (MEST). 
The work of S. Sharpe is supported in part by the US 
DOE grant no. DE-FG02-96ER40956. J. Bailey, FL- 
J. Kim and W. Lee thank the University of Washington 
for hospitality while this work was begun. 



[1] 



J. Laiho, E. Lunghi, and R. S. Van de Water, Phys.Rev., 
D81, 034503 (2010)| |arXiv:0910.2928 [hep-ph]| ~ 

S. Durr, A. Juttner, L. Le llouch, 
C71, 1695 (2011)| 



[2] G. Colangelo, S~D urr : 

H. Leutwyler, et al, Eur.Phys.J. 
|arXiv: 101 1.4408 [hep-latj | 
[3] M. Ciuchini, V. Lubicz, L. Conti, A. Vladikas, A. Donini 
et al, JHEP, 9810, 008 (1998), erratum added online 

Mar/29/2000, |arXiv:h ep-ph/9808328 [he p-ph] | 

[4] A. J. Buras, S! Jager, and J. Urban, |Nucl.Phys. , B605 



600 (2001)| |arXiv:hep-ph7 0102316 [ hep-phj | 
[5] M. Ciuchini, E. Franco, D. Guadagnoli, V. Lubicz 
V. Porretti, et al, | JHEP, 0609, 013 (2006)| |arXiv:hep"- 
ph/0606197 [hep-ph] 



arXiv:1201.1302 



6| A. J. Bu ras and J. Girrbach, (2012), 
[hep-ph]] 

M. Bona et al (UTfit Col laboration), JHEP, 0603, 
(2006)| |a rXiv:hep-ph/0509219 [hep-ph]| 
Isidori, 



[7] 



080 



[8j G. Isidori, Y. Nir, and G. Perez, |Ann7Rev.Nucl.Part.Sci 

60, 355 (2010)1 |arXiv:1002.b 900 [hep -phj | 
[9j A. Hasenfratz and F. Knechtli,|Phys. Rev., D64, 034504 
(2001)1 |arXiv:hep-lat/0103029 [hep-latj | " 
[10J C. W. Bernard, T. B urch, K. Orginos, D. Toussaint, 
T. A. DeGrand, et a/.,|Phys. Rev., D64, 054506 (2001)| 
|arXiv:hep-lat/0104002 [hep-latj| 
[11] W.-J. Lee and S. R. Sharpe, |Phys.Rev ., D60, 114503 



(1999), arXiv:hep-lat/9905023 [hep-lat 



12| C. Bernard fMILC Colla boration), Phys.Rev., D65 
054031 (2002)| [arXrv:hep-lat/0111051 [hep-lat~ ~ 



M 

14 



[18 



w 

[20 

r 



[21 



22 



D83, 



C. Aubin and C. B ernard, Phys.Rev ., D68, 034014 

(2003) [ |arXiv:hep-lat/0304014 [hep-latj| 
J. Kim, W. Lee, and S. R. Sharpe, Phys.Rev 
094503 ( 2MT)|[arXiv:1102.1774 [hep-latj| 
R. S. Van de Water and S. R. Sharpe, |Phys. Rev., D73, 
014003 (2006)] |arXiv:hep-lat7 0507012 

D. Becirevic and G. Villadoro, |Phys.Rev., D7 0, 094036 

(2004) , ar Xiv:hep-lat /040 8029 [he p-latj| 
F. Gabbiani, E. Gabrielli, A. Masiero, and L. Silvestrini, 
|Nucl.Phys., B477, 321 (1996)| |arXiv:hep-ph/9604387| 



[hep-ph] 

C. Allton, L. Conti, A. Donini, V. Gimenez, 
L. Giusti, et aZ.,|Phys.Lett., B453, 30 (1999)||arXiv:hep-| 
lat/9806016 [hep-latj [ " 
al. 



T. Bae et al, |Phys. Rev, D82, 
|arXiv:arXiv: 1008.5179 [hep-latj | 



114509 



(2010) 



and G. Martinelli, 
|arXiv:hep-lat/9910017 



A. Donini, V. Gimenez, 
Phys.Lett ^ B470, 23T 
[hep-lat] 

G. Kil cup and S. R. Sharpe, |Nucl.Phys., B283, 493] 
(1987)| 



G. Kilcup, R. Gupta, and S. R. Sharpe, Phys.Rev., D57, 



1654 (1998), arXiv :hep-lat/9707006 [hep-latj| 
M. Ciuchini, E. Franco, V. Lubicz, G. Martinelli, 



[24] 



I. Scimemi, et al, 


Nucl.Phys., B523, 501 (1998) 


arXiv:hep-ph/9711402 [hep-ph] 





397 (2000) , arXiv:hep-ph/0005183 [hep-ph] 



16 



[25] R. Gupta, T . Bhattacharya, and S. R. Sharpe, Phys. 

Rev., D 55, 4036 (1997))|arXiv:hep-lat/9611023 [hep-lat] | 
[26] H. Kluberg-Stern, A. Morel, O. Na poly, and B. Peters- 

son, |Nucl.Phys., B220, 447 (1983)| 

[27] D. Verstegen, |Nucl.Phys., B249, 685 (198571 

[28] W. Lee and Sharpe, |Phys. Rev., D68, 054510 

I (2003)| |arXiv:hep-iat /0306016. 

[29] J. Kim, W. Lee, and S. R. Sharpe, (private communi- 
cation). 



[30] S. R. Sharpe and R. S. Van de Water, 


Phys.Rev., D71, 


114505 (2005) arXiv:hep-lat/0409018 


hep-lat] 



[31] C. Bernard, |Phys.Rev., D73, 114503 (2006)] |arXiv:hep- 

lat/0603011 [hep-lat], 
[32] C. Bernard, M. Golterman, and Y. Shamir, Phys. Rev.,| 

D77, 074505 (2008)| |arXiv:arXiv:0712.25 60 [hep-latjp~ 
[33] S. R. Sharpe and N. Shoresh, [ Phys.Rev., D62, 094503[ 

(2000)[|arXiv:hep-lat/0006017 hep-lat] | " 
|34| S. R. Sharpe and N. Shoresh, Phys.Rev., D64, 114510| 

(2001) , ar Xiv:hep-lat/010 8003 [he p-lat] | 

[35] C. W. Bern ard, T. Draper, A. Soni, H. P olitzer, and 

M. B. Wis e |Phys.Rev., D32, 2343 (1985)| 
[36] A. Roessl, |Nucl. Phys., B555, 507 (1999)] |arXiv:hep- 

ph/9904230 

[37] C. Alrton~! et al. (RBC-UKQCD), |Phys. Rev., D78, 
114509 (2008)[[arXiv:0804.0473 [hep-lat] | 



